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ABSTRACT 


The unknown boundary conditions at the machine tool 
bed-carriage interface make it difficult to analyse the joint. At 
the contact nodes both the displacements as well as forces are 
unknown. Depending on the contact forces and the coefficient of 
friction, the nodes will either slip or adhere to each other 
during deformation. The deformation may also render some nodes 
out of contact. 

In this work, a solution procedure to analyse the contact 
problems with varying coefficient of friction is developed using 
the finite element technique. This is further used to analyse the 
V-slide joint of a planer bed for studying the deflections and 
stresses for different values of coefficient of friction, and 
different points of application of the applied force. 



CHAPTER I 


INTRODUCTION 

1.1 MACHINE TOOL JOINTS: 

Machine tool joints may be classified as fixed joint such as 
bolted connections, shrink fit assembles, etc. or sliding joints 
such as various types of guideways. Various forces including the 
frictional forces are transmitted across the interface which cause 
elastic displacements. The interface displacements influence the 
overall displacements of the machine tool members. Tilting of the 
slides due to deformation affects the operating properties of the 
si ideways . 

It is essential to know the deflections as they have a 
profound influence on the machine tool operational accuracy. From 
strength point of view, it is necessary to know whether the 
stresses are below the safe levels. In particular shear stress 
plays an important role in wear of contact surfaces. 

1.2 LITERATURE SURVEY 

1.2.1 Methods f or Determining Contact Deformations: 

For determining the pressure distribution at the contact 
surface and the deformations due to external load N. Back, M. 
Burdekin and A. Cowley [1] proposed an iterative approach. In this 
approach, the compliance of the contact surface is assumed 
nonlinear as given by the relation: 
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X 


( 1 . 1 ) 


= C (100 P) ro 

where 

^ = approach of surface in ^m 

® ~ constants depending on material and surface finish 
characteristics 
p = pressure 

The contact between the components is simulated using springs or 
plates with properties which are taken as functions of surface 
roughness, pair of materials and pressure distribution. . This 
method however has limited application as it relies on the 
empirical relation (1.1). Further, it does not take into account 
the effect of friction. 

Another approach of calculating stresses and deformation in 
machine tool joints is to analyse the joint as a contact problem. 
Various classical solutions for this special class of solid 
mechanics problems are available. In three dimensions, the 

problem of contact of two elastic bodies was first formulated and 
solved by Hertz [2] under several restrictive assumptions. The 
Hertz problem was reduced to an integral equation and solved by 
Mushkhel ishv i 1 i [3] and Gladwell [4]. However, the complex 
geometries of machine tool joints are extremely difficult to model 
in the conventional way. Here finite element methods plays an 
important role in analysing such complicated geometries and 
boundary conditions. 

1 . £s. 2 Finite Element Methods for Contact Problems: 

The earliest application of finite element method to contact 
problem appears to be in the paper of Chan and Tuba [5]. S. Ohte 
[6] proposed a method in which the contact interface is classified 
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into adhering and sliding states according to reaction force 
acting on the contact surface and the friction coefficient. 
Francavilla and Zi enki evicz [7] solved the frictionless case of 
contact problems as a quasilinear problem. Hung and Sauxe [8] 
used finite element method and mathematical programming techniques 
to solve frictionless contact problems. Here appropriate 
linearization of contact conditions is adopted and solved as an 
non-iterative optimisation problem. Sachdeva, Ramakrishnan and 
Natrajan extended the method of Francavilla and Zienklewicz to 
take into account the force bot.idary conditions [9] and frictional 
effects under proportionate loading [10], Torstenfelt [11] 
proposed a general purpose finite element computer program to 
solve contact problems with friction. P.V. Kishore [12] proposed 
the minimization of dissipation energy approach for solving 
contact problems. Here the extent of slip is determined by the 
minimisation of the dissipation energy function. 

The approach suggested by Sachdeva and Ramakrishnan [10] has 
been modified in the present work. They have used flexibility 
matrix, whereas in the present approach the contact conditions are 
applied to the stiffness matrix. Besides this, the matrix 
substructuring technique proposed by P.V. Kishore [12], is 
utilized for optimum storage and computational efficiency. 

1.3 OBJECTIVE AMD SCOPE OF PRESENT WORK: 

The present work describes an iterative numerical procedure 
to solve two dimensional elastic contact problems with friction 
using the conventional stiffness matrix. This stiffness matrix is 
condensed using a substructuring algorithm to eliminate all the 
nodes except those which are likely to be in contact and those 
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where there are external forces acting. To begin with a set of 
nodes is assumed to be in contact. Newton’s third law is applied 
to nodal forces and the compatibility conditions are applied to 
both the normal and tangential nodal displacements. In the 

subsequent iteration, the contact conditions are dropped for the 
nodes going out of contact and for the slipping nodes, the 
displacement compatibility is applied only in the normal direction 
while the slip condition is applied in the tangential direction. 
Stresses at the contact surface are obtained by averaging the 
stresses over tne bodies in contact. The iterative scheme 
converges very fast. 

The method presented here is first applied to a simple 
problem with known solution to test its applicability, followed by 
the planer tool bed-carriage V-joint problem. The two bodies in 
contact are the portions of carriage and the bed. Although the 
cutting force has both the components, only the vertical component 
is considered in analysis. It is assumed that the horizontal 
force will not affect much the pattern of deformation/stresses in 
the contact region. The above assumption helps us in idealising 
the configuration as a plane strain problem whereby only a thin 
slice of the body needs to be analysed. Only half the 
carriage-bed structure is considered by assuming symmetric loading 
conditions and imposing suitable boundary conditions at the line 
of symmetry. The three cases of coefficient of friction account 
for the dry metal to metal contact and boundary lubrication 
conditions . 

The limitations of the method may be summarized as follows: 
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(a) The contact conditions for the two bodies are applied only at 
the nodes, whereas ideally they have to be applied over the 
entire contact length. 

(b) The cutting forces have been assumed to act at a point 

whereas in actual conditions there is a finite width of cut. 

(c) Body forces have been neglected in the present analysis as 

their magnitude is small as compared to the cutting force. 

(d) Incremental loading has not been considered. But in starting 
operations, the magnitude of cutting force actually increases 
from zero to a peak value. 

Since the problem has been idealised as a 2-D problem, the 

results are valid only in the middle portion of the carriage where 

plane strain conditions are more or less valid. However, the 

actual problem is 3-Dimensional. The effect of the horizontal 
forces can also be studied only in a 3-Dimensional analysis. 

1.4 PLAN OF THE THESIS: 

Chapter II deals with the finite element formulation of a 
contact problem. It describes the details of the substructuring 
and the application of the contact conditions. The results and 
discussion are given in Chapter III. The work concludes in 
Chapter IV with the suggestions for future work. 
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CHAPTER II 


FINITE ELEMENT ALGORITHM 


This chapter 

presents 

the 

f inite 

element formulation 

o£ 

elastic bodies in 

contact 

£ or 

plane 

sTain conditions. 

In 


particular, it describes the matrix substructuring technique where 
a condensed stiffness matrix is obtained from the stiffness 
matrices of bodies in contact in skyline form. Besides, the 

contact conditions algorithm is also discussed. 

2. 1 FEM FORMULATION FOR PLANE STRAIN 

Uhen a cylindrical body is subjected to support conditions 

and external forces which have no components along the axis of the 

cylinder and the remaining components do not change along the 

axis, the body is said to be in a state of plane strain. If the 

2-coordinate is taken along the axis of the cylinder, then the 

displacement, stress and strain components (with the exception of 

<y 1 associated with the 2-direction are zero and the remaining 

z 

components are independent of 2. Then it is sufficient to analyse 
only a thin slice of the body. 

Figure 2.1 shows a slice of thickness h of a cylindrical body 
in a state of plane strain. Let ft be the region occupied by the 
cross-section of the slice and r be its boundary. The part of 
the boundary on which the external forces act is denoted by T q . 
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The relation between the non-zero strain and displacement 
components is given by 

3\i 3v Sxi 3 v r o i > 

45 = - 5 — ; 45 = ; 45 = — - + — - t Z . 1 ; 

X 3x f y <*y’ xy 3y 3x 

If we write the strain and displacement components as 


{45} 


“{«*). w ■{:} 


( 2 . 2 ) 


xy 


then this relation becomes 


{ 45 } 


3 

3x 


3 

Sy 


0 

3 

3 

3x 


{“} 


(2.3) 


If the stress components & , a , r are also arranged in a vector 

x y xy 


f orn 


<*■> 


•{>) 


(2.4) 


xy 


Then the stress-strain relation for a linearly elastic material 
under plane strain conditions, becomes 
{<*> = [D] («} 

where [D] Is given by 


(2.5) 


E (1 - **) 

CD3 = (i+v) TT ^fe) 


w / ( 1 -V ) 
0 


v/(l-v) 0 

1 0 

0 (l-2v)/2(l-iO 


( 2 . 6 ) 
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where E is the Young’s Modulus and 


v is the Poisson’s ratio of the material. 
The strain energy is given by the expression 

U = J j {«} T {&) h dA 


Cl 


= J j <«) T < D> <«} h dA 


Cl 


Let [f > 
s 


■{':.} ■ 


{B f } 


■C:) 


(2.7) 


( 2 . 8 ) 


represent the applied traction on and the body force (per unit 
volume) respectively, then the work done by the external forces is 
given by 

Up = J \ (B) T {u} h dA + J {f s } T (u> h dA 

ci r 

q 

(2.9) 

The expression for potential energy of the system is given as 


rr = U - U 
P P 


Hence 


= j J <<=> T (D) {<=} h dA - J \ (B f > T {u> h dA 


Cl 


Cl 


- J {f } T <u} h dA 


( 2 . 10 ) 


Let the domain be divided into N g number of elements (Figure 
2.2). Over a typical element Q g , let 
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<U> = 


N 1 0 
0 N„ 


N 0 
n 

0 N 

n 


u. 


n 


n 


[N] <u} 


nodal 


( 2 . 11 ) 

be the approximation for the displacement vector. Here u® , ... 

e & 

u^, v n are the values of the displacements at the n nodes of the 
element while ...N^ are polynomial functions of the coordinates 
x and y, called the shape functions. From equation (2.3) we get 


{«} = [B] (u> 


nodal 


( 2 . 12 ) 


where the matrix [B] contains the derivatives of the shape 
functions. Substituting equations (2.11) and (2.12) into (2.10), 
the expression for potential energy can be written as 


n 

P 


N T 

L 7 I '“’Jodal I B]T 1D! tB1 '“’nodal hdA 
e_1 a 

e 

N T 

- 2" h J '“ ’nodal "*’ T < B f> dA 

‘- 1 a 

e 

N b T 

- I»J '“’nodal [N ' ]T <B a> dl 

r b 


(2.13) 


Note that when the domain is divided into area elements, the 
boundary F automatically gets divided into line elements. N fe is 
the number of such line elements, r fc is the domain of a typical 
line element b and 
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is the approximation for the displacement vector over it. Setting 
the first variation of n to zero, we get 

p ^ 

t K] {U} = { F } 

where 

N 

e 

[K] = ^ [k] 6 (Stiffness Matrix) 


[k] 6 = J [B] T [D] [B] h dA 

e 

N N k 

e d , 

{F> = ^ {f} 6 + ^ -{ f ^ (Force Vector) 

e “ 1 b-1 

{k} 6 = J [ N] T { B f } h dA e 

ft 

e 

{f > b = J [N] >T <f > h dA 

J S 6 


(f) e = J tN] T <B f > h dA e (2.15) 

r 

e 

and { U > is the vector of nodal displacements for the entire 

domain. In eq. (2.14) it is understood that the element matrix 

a eh 

[k] and element vectors {f} and {f> are first expanded to the 

full size before adding. 
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2.2 FORMULATION FOR CONTACT PROBLEMS: 


Two bodies A and B which are in contact with each other are 
shown in Fig. 2.3. The finite element equations for these bodies 
are 


[K a ] <u a ) = (F a > 

[K b ] <U b > = <F b > (2.16) 

where [1^], [Kg], {U A >, {U B >, { F^} , {Fg} are the stiffness 

matrices, displacement vectors and force vectors of bodies A and B 
respectively. The displacement and force vectors include the 
contact displacements and forces also. These equations can be 
combined into a single equation: 


[K a ] 


[K b ] 


f< U A>) 




f<F A > 


1 ( f b } 


or [K] {U} = {F} (2.17) 

The equation (2.17) cannot be solved in the given form as at 
the contact nodes, neither the displacements nor the forces are 
known. To overcome this difficulty contact conditions are used. 


2.3 CONTACT CONDITIONS: 

There are two set of contact conditions, one on the unknown 
contact forces and the other on the displacements. 

Let the i-th node ofbody A be in contact with the j-th node 
of body B (Fig. 2.3). Let ’n’ and 's’ denote the normal and 
tangential directions respectively at the contact node. 
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F 



Fig. 2-3 Two bodies A and B in contact. 



Fig.2-4 Re-arranged matrix. 
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The force conditions given by Newton’s third lav are: 


< f a> 

< f a> 


in 

is 


< F B> 

< f b> 


jn 

js 


The displacement compatibility conditions are 


(2.18) 

(2.19) 


(U. ) . = (U n ) . + « 

Am B'jn 


(Vi. 


<Vjs 


( 2 . 20 ) 

( 2 . 21 ) 


where « denotes the initial clearance between the nodes which come 
in contact . 

However, the tangential displacement compatibility (2.21) is 
valid only when the nodes do not slip. Slip occurs when the shear 
force exceeds or equals the frictional resistive force. 
Frictional resistive force is dependent on the coefficient of 
friction and the normal force at the node. Therefore, if slip 

occurs equation (2.21) is replaced by 

l (F A>j=l ’ " l (i Vjnl (2 - 22:i 

This condition can be applied to body B also. 

In order to apply these conditions to the system of equations 
(2.17), they have to be expressed in terms* of the ’ x ’ and ’y’ 
components. Resolving ’n* and ’s’ components along the x and y 
axes, we get (See Fig. 2.3). 


'Vin * '‘Vi. £in “ + (F A>iy C °° a 

(F A 5 is = +(F A ) x* Cos “ + ‘Viy Sin “ ( ' 2 ' Z3:> 

where is the angle which ' s ’ -direct ion makes with x-axis. 

Similarly (F fi ) jn , (F B ) jg , (U A >in* C^is' (Vjn and (U B )js 
can be expressed in terms of x and y components. Substituting eq. 
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(2.23) and similar relations in equations (2.18) - (2.22). Ue then 
eet , 

~ CF A 5 ix Sin °* + (F A 5 iy Cos “ " ( F B >j x Sin « + C F B ) jy Cos « = 0 

(2.24) 

(F A ) ix Cos « + (F A ). y Sin « + (F fi ). x Cos « + (F B ) jy Sin <* = 0 

(2.25) 

Sin “ " CViy Cos « - CFg) jx Sin ,« + (F fi ) jy Cos « = 0 

(2.26) 

ix Cos 01 + ^Viy Sin « - CU R ) jx Cos « - (U fi ). y Sin a = 0 

(2.27) 

((ss Cos cn + fj sn Sin cx) (F n ) . + (ss Sin a - #j sn Cos «)(F B ) . = 0 

c jx a jy 

(2.28) 

where ss is the sign of (F n ). and sn that of (F n ). . Note that 

o j s d j n 

Eq. (2.27) is applied for no slip condition while eq. (2.25) is to 
be used when there is slip. 

2.4 Application of Contact Conditions: 

It can be seen that corresponding to every contact pair (i, 
j), there are four equation in the set (2.17). Two equations 
(21—1 ) - th Eq. for x-directlon and (2i)-th for y-direction for the 
body A and the remaining two, (2j-l)-th Eq. and (2j)-th, for the 
body B. As stated earlier, neither the displacement associated 
with these nodes nor the forces (i.e. the right hand side of these 
equations) are known. Therefore, these equations cannot be solved 
as such. However, for every contact pair, there are four contact 
conditions. The strategy is to replace these equations with 
unknown right hand sides by the contact conditions which have the 
known right hand side. The matrix operations required for this 
process are described below. 
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First the (21 1 ) - th Eq. of body A is replaced by Eq. (2.24). 

Since and (F^)^ are equal to the left hand sides of 

(2i-l)-th and (2i)-th Eq . of body A and (F D ). and (F_) . are 

d jx 15 j y 

equal to the left hand side of (2j-l)-th and (2j)-th Eq. of body 
B, the condition states that 

(- Sin ot) ( 2 i - 1 )th Eq . of body A 

+ (Cos ot) (2i)-th Eq. of body A 
+ (-Sin oi) (2j-l)-th Eq . of body B 
+ (Cos «) ( 2 j ) - th Eq. of body B = 0 

Thus after replacement, the (2i-l)-th row of the [K^] becomes 

equal to 

(- Sin oi) (2i - l)-th row of [K^] 

+ (Cos oi) ( 2 i )-th row of [ ] 

+ (-Sin oi) (2 j-l)-th row of [K fi ] 

+ (Cos «) (2j)-th row of [Kg] 

The (2i-l)-th row of [ ] is set to zero. 

Next, the (2i)-th Eq. of body A is replaced by the Eq . 

(2.25). The matrix operations for this replacement are similar. 
Now, the (2i)-th row of [K^] is replaced by the combination of 
(2i-l)-th and (2i)-th rows of [ K^] and (2j-l)-th and (2j)-th rows 
of [K n ], the coefficients of the combination given by the Eq . 

(2.25) . The (2i)-th row of [F A ], becomes zero. Note that the Eq . 

(2.25) involves the original (2i-l)-th row of [K^] and not the 

modified one. Thus the (2i-l)-th row of [ ] has to be stored 

before applying the first contact condition. 

The third condition is applied in the end, the reason for 
which is explained later. The fourth contact condition, discussed 
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hereby. 

is used to replace the (2j-l) 

-th Eq . 

of 

body B. Eq . (2 

.27) 

is used 

for the no slip condition. To 

apply 

this condition. 

the 

(2j-l)- 

th row of [Kg] is made zero 

except 

f or 

the element 

in 

C 2 i-1 ) , 

( 2 i ) , ( 2 j — 1 ) and ( 2 j ) columns 

whi ch 

are 

respectively 

Cos 

oi. Sin 

<*, (-Cos «) and (-Sin ot) . The 

( 2 j — 1 ) 

-th 

row of {Fg> 

is 


made zero. Thus the application of displacement contact condition 
is more simpler than that of force constant condition which 
involves certain row operations. In case of slip, Eq . (2.28) is 

used. This condition is similar to the first two conditions 
except that it involves c ily two equations, i.e., (2j-l)-th and 

(2j)-th equations of body B. Here, the (2j-l)-th row of [Kg] is 
replaced by the linear combination of (2j-l)-th and (2j)-th rows 
of [Kg], the coefficients of the combination being given by Eq . 
(2.28). the (2j-l)-th row of [Fg] is then made zero. 

Finally, while applying the third contact condition the 
( 2 j ) -th row of [Kg] is made zero except the elements in the 
columns (2i-l), (2i), ( 2 j — 1 ) and ( 2 j ) which are given by (Sin ot) , 
(-Cos oi) , (-Sin «) and (cos «) respectively. The (2j)-th row of 
[F_] is replaced with clearance («). 

JD 

In case the third contact condition is applied before the 
fourth, (2j)-th Eq. of body B gets modified. However, for slip, 
the original (2j)-th Eq . is required for the fourth condition. So 
by applying fourth earlier, we avoid storing (2j)-th row of [Kg]. 

Note that application of the contact conditions makes the 
coefficient matrix unsymmetric. 

2. 5 SUB -STRUCTURING* 

An iterative scheme involving large unsymmetric matrices 
requires a large storage and a considerable amount of computing 
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time. To reduce both storage and computing time, the stiffness 
matrix is stored in skyline form and then decomposed to a much 
smaller size using sub-structuring. 

The first step in substructuring is to rearrange the rows and 
columns of the global stiffness matrix as shown below: 


I, I 

K I , II 

K I , III 


r U I ' 


f 

F I ^ 

II, I 

K II , II 

K II, III 


U II 

y = « 


F II > 

III , I 

K I 1 1 , 1 1 

K III , III 


‘ U III - 


L 

F I I I - 


(2.29) 

Here the suffix I denotes the degrees of freedom associated with 
the contact nodes and the nodes at which non-zero forces are 
specified, II denotes the degree of freedom associated with the 
free surface and the internal nodes, while III denotes the degree 
of freedom associated with the nodes at which displacements are 
specified. If the prescribed displacements are zero then 


F n - 0 


Then the equation (2.29) reduces to 


K I,I U I + K I,II U II = F I 
K II,I U I + K I I , I I U II = ° 

Eliminating U n from the above equations we get 

v = F 

i,I Uj I 


(2.30) 


(2.31) 
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where 


•i,i 


(K i , i " K i 


n K n , ii K n 


To obtain the matrix j , first Kj j is extracted 
form from the global stiffness matrices of bodies A 
skyline form. Next Kjj j j. is obtained in skyline form 
global stiffness matrices. Then (K”* n Kj j) is 
column-wise. If { C } ^ is its i-th column, then it can be 
solving the system 


(2.32) 

in full 
and B in 
from the 
obtained 
found by 


II, 11^ ^II,lh-th column 

Thus the inverse of jj is not explicitly found. Then the 
elements of the matrix Kj jj Kjj jj Kjj j are obtained by 
multiplying the rows of j j with the appropriate column of 
Kjj ij Kjj j. In the last two steps, the rows of Kj jj or the 
columns of Kjj j. are generated from the skyline form as and when 
required. Finally ~ is obtained using Eq . (2.32). 

Substituting this leads to a very small set of equations 
(2.31) which is solved after applying the contact conditions. 
Displacements Ujj are obtained by back substituting Uj in the 
equation 


U 


II 


~ K II,II K II,I 


(2.33) 


Once all the displacements are known Fj can be readily found using 
the third set of Eq. (2.29). 


2. 6 Program Algorithm 

The set of equations (2.31) are solved iteratively by 
applying the contact conditions in the following way. 

1 A set 0 f nodal points is assumed to be in contact (it should 


Cr\r .,* K L r;pn W 
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include all possible contact points) and equations (2.31) are 
solved after applying the contact conditions for no slip. 

2. The sign of the normal force at each node in the contact one 
is checked. The normal nodal forces cannot be tensile. Hence all 
those nodes where the normal force turns out to be tensile are 
deleted from the possible contact zone in the next iteration. 

3. The ratio of shear force to normal force, at each node in the 
contact zone is checked. Slipping occurs when the nodal shear- 
force exceeds or equals the product of nodal normal force and the 
coefficient of friction. For the nodes which il ip , the condition 
(2.26) is replaced by eq. (2.27) in the next iteration. 

4. Equation (2.31) are solved by repeating steps 2 and 3 until 
all the normal nodal forces in the contact zone come out to be 
compressive and the ratio of the shear force to the normal force 
for all the nodes in the contact zone becomes either equal to or 
less than the coefficient of friction. 


2.7 EVALUATION OF STRESSES; 

To evaluate the stress components at a point, first the nodal 
displacement vector ° f the elefflent to which the point 
belongs is extracted from the solution vector (U). Stresses are 
then calculated using the Eqs . (2.5) and 2.12) i.e. 

<«.} = [D] <B> <«>' odal (2-34) 

Here we need stresses in the contact region which is the 
common boundary between bodies A and B. Normally the 
approximation chosen for the displacement vector Ceq- 2.11) is 
such that the stresses do not become continuous across the inter 
element boundaries. Thus the stress components at a point on the 
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contact surface will have different values when it is considered 


as a point belonging to body A first and then to body later. 

Therefore we adopt the following procedure on the contact surface. 

First we choose a set of points which are different from the 

contact nodes. The stress components <y , o* and r at these 

x y xy 

points are determined using the relation (2.34) by assuming that 
they belong to body A first and then to body B later. Next the 
X~Y components of the stress vector on the contact surface are 
calculated using the relation 



where n and n on the components of the unit vector along normal 
x y 

direction n with respect to xy axes. Then, normal and shear 

component of the stress vector are given by 


f = f n + f n 
n xx y y 


f = f s + f s 
s xx y y 


( 2 . 36 ) 


where s and s are the components of the unit vector along the 
x y 

tangential directions with respect to xy axes. 

Finally, the magnitudes of the normal and shear components of 

the stress vector (f } are averaged over the bodies A and B to 

s 

give the contact stresses. 

It may be noted that the maximum tensile stress for the 
planer bed-carriage problem is directly determined by using the 
relation (2.34), where <u)* odal ia derived for the bottom centre 
element of the carriage from the solution vector {U>. 


21 



CHAPTER III 


RESULTS AND DISCUSSIONS 

"this chapter, first the problem of flat plate over flat 
plate for which the solution is known is solved to verify the 
applicability of the program. Then the problem of V-slide of a 
planer bed-carriage is considered. All the materials are assumed 
to be isotropic and linearly elastic. Body forces are neglected 
and axial thickness ’h' is taken as unity. In both the problems, 
four noded isoparametric element is used. The program developed 
is in general for unsymmetric problems, however the 
simplifications arising out of symmetry have been considered. 

3.1 FLAT PLATE OVER A FLAT PLATE 

Contact between two cast iron plates of same material with a 
point load at the mid-length on the top plate is analysed here. 
The configuration is shown in Figure 3.1(a). 

For the finite element mesh, body 'A’ is divided into 60 
elements with 84 nodes and body ’B’ is divided into 96 elements 
with 125 nodes. One half of this is shown in Figure 3.1(b). 

For computation the following data is used: 


Force at the middle. 

P = 

1960 N/mm 


Young’s Modulus, 

E = 

0.931 x 10 5 

N/mm^ 

Poisson’s ratio. 

v = 

0 . 25 


Coefficient of friction 

V = 

0.002, 0.20 

o 

O 

*k 
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Fig. 3-1 (a) Flat plate over a flat plate configuration. 
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Fig. 3 . 1(b) One half of the finite element mesh and 
deformed shape. 
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Normal Stress, N/mm 



Semi-contact length, mm 

Fig. 3- 1(c) Variation of normal stress along the contact 
length for various values of p. 



To begin with 21 nodes are initially assumed to be in 
contact. In the converged solution only 13 nodes are found to be 
in contact for p = 0.002 and (j = 0.20. For ^ = 0.40 the contact 
length extended to 15 nodes. 

The deformed shape of the two bodies in contact for (j — 0.20 
is shown to an enlarged scale in Figure 3.1(b). The normal 
pressure distribution for various values of is shown in Fig. 
3.1(c). It is found that the contact zone increases with increase 
in the value of the coefficient of friction and the maximum normal 
pressure decreases with increase in friction. However, the 
magnitude of this reduction is not very prominent. This may be 
due to the fact that the point of maximum pressure is under the 
direct influence of the applied load. 

The results obtained are found to conform very closely with 
that of Ramesh [13]. The results of Sachdeva and Ramakrishnan 
[12] were validated by Ramesh [IS] who had solved the same 
probl em . 

3.2 PLANAR BED - CARRIAGE CONTACT PROBLEM 

A planer has a moving carriage which slides in V-shaped 
slidway of the bed. The cross-section of the carriage-bed 
configuration Is shown in Figure 3.2(a). The bed and carriage 
dimensions were obtained from the machine no. 31, Central 
Workshop , I IT, Kanpur. The vertical cutting force data has been 
used from the reference [14] which takes into account the heavy 
and rough machining conditions encountered in planing. The force 
location has been varied from point (1) to (4) (Refer Figure 3.2 
(b)). Since the structure is symmetric in geometry and loading 
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f ig3 2 (r) V- joint of .1 Planer. 
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fig 3.2(b) Finite element me* 


conditions, only half the cross-section is analysed. The symmetry 
conditions have been accounted by imposing u-displacements and 
shear stress equal to zero in the centre plane. Each case has 

been analysed for three values of friction accounting for the dry 
and boundary lubrication conditions. 

The joint and the maximum structural deflections, contact 
length, the normal and shear contact stresses and the maximum 
tensile stress in the structure, have been determined to study how 
far the strength and deflection criteria of safe design have been 
incorporat ed . 

The finite element mesh is shown in figure 3.2(b). Carriage 
has been divided into 103 elements with 145 nodes while the bed 

has been divided into 56 elements with 75 nodes. Since the bed 

and carriage are made of cast iron, the following properties are 

us ed : 

Vertical force, P = 300 N/mm 

Young’s Modulus, E = 0.931 x 10^ N/mm^ 

Poisson's ratio, v = 0.25 

Coefficient of friction, v = 0.002, 0.20, 0.40 

2 

Strength in compression, ^cu = 500N/mm 

2 

Strength in tension, <^ tu = 165 N/mm 

2 

Strength in shear, o- su = 50 N/mm 

The general deformation pattern of the structure is shown in 
Figure 3.3. It corresponds to the force location at point (4). 
Variation of the normal deflections along the contact length for 
various values of fj and force locations is shown in Figure 3.4. 
The corresponding tangential deflections for the bed and the 
carriage are shown in figures (3.5) and (3.6) respectively. The 
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f ig.3 7(a) Variation of normal and shear stresses along the 

contact length for various values of y 
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Node no: 



r— 1 1 i J 0 

1 2 345676 9 10 



fig. 3.7(b) Variation of normal and shear stresses along the 
contact length for various values of jj 



FORCE AT 3 


Node no : 



N/4nm^ 



contact length for various values of JJ 



f ig, 3.7(d) Variation o f norma! and sheer stresses aiong 
Contact length for various values of |J 



variation of the magnitude of normal and shear stresses along the 
contact length for various values of p are shown in figures 
3.7(a), 3.7(b), 3.7(c) and 3.7(d) for four different locations of 
the applied force. For the carriage shear stresses on both the 
parts AB and CD act downwards. Then for the bed, they will be in 
the upwards direction. 

The observations for the four force locations are as follows: 
Caae 1 : Force at location ( 1 ) . 

All the nodes slip initially but remain in contact till 
convergence. The solution converges in three iterations. The 
bottom centre node of the carriage deflects down vertically by 
0.72 /urn. The tensile stress in the neighbourhood of this point is 
-2.3 N/rnm 2 . 

The salient data for the contact zone is given in Table A. 


Table A: Case (1) Data 


Value of 

P 

Node 

No. 

Max . 
Normal 
def . 

xl0~ 2 

(mm) 

Node 

No. 

Max . 
Tang, 
def . 

xlO -2 

(mm) 

Node 

No. 

Max . 
Coropr . 
stress 

(N/flim^ ) 

Node 

No. 

Max . 

Shear 

stress 

(N/mm^ ) 

0.02 

10 

0 . 435 

10 

0.384 

6 

12.0 

6 

1.6 

0.2 

10 

0.431 

10 

0.382 

6 

11.0 

6 

1.3 

0.4 

10 

0.431 

10 

0.382 

6 

8.2 

6 

1.0 


Caae 2 : Force at location (2) 

All the nodes slip initially but remain in contact till 
convergence. The solution converges in three iterations. The 
bottom centre node of the carriage deflects down vertically by 
0.83 pm. The tensile stress in the neighbourhood of this point is 
-0.7 N/mm 2 . 


38 



The salient data for the contact zone is given in Table B. 


Table B: Case (2) Data 


Value of 

Node 

No. 

Max . 
Normal 
def . 

xl0~ 2 

(mm) 

Node 

No. 

Max . 
Tang, 
def . 

xl0~ 2 

(mm) 

Node 

No. 

Max . 
Compr . 
stress 

(N/mm 2 ) 


Node Max. 

No. Shear 

stress 

( N / mm ^ ) 

0.02 

10 

0.436 

10 

0.384 

6 

11 . O' 


6 

1 . 4 

0.2 

10 

0.435 

10 

0 . 383 

6 

9.8 


6 

1 . 2 

0.4 

10 

0.435 

10 

0 . 381 

6 

8.6 


6 

1 . 1 

Case 3: 

Force 

at location (3) 

1 . 






All 

the nodes slip 

initially and node 

1 goes 

out 

of contact 

before convergence. The 

solution converges in 

four 

iterations . 

The bottom centre node 

of the 

carriage deflects down 

i vertically by 

1.55 /urn. 

the tensile stress in the 

neighbourhood 

of 

this point is 

31.3 N/mm 2 . 









The 

salient data for the 

contact zone 

is given 

in 

Table C. 



Table C: 

Case 

(3) Data. 




Value of 

Node 

No. 

Max . 
Normal 
def. 

xl 0~ 2 
(mm) 

Node 

No. 

Max . 
Tang, 
def . 

xlO -2 

(mm) 

Node 

No. 

Max . 
Compr . 
stress 

(N/mm 2 ) 


Node Max. 

No. Shear 

stress 

(N/mm 2 ) 

0.02 

10 

0.440 

10 

0.385 

5 

11.0 


5 

1.4 

0.2 

10 

0.439 

10 

0.384 

5 

9 . 6 


5 

1.3 

0.4 

10 

0.439 

10 

0.383 

5 

8.3 


5 

1.1 


Case 4^ Force at location (4). 

All the nodes slip initially and nodes 1 and 2 go out of 
contact before convergence. The solution converges in four 
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iterations. The bottom centre node of the carriage deflects down 
vertically by 1.59 pm. the tensile stress in the neighbourhood of 
this point ia 48.7 N/mm 2 . 

The salient data for the contact zone is given in Table D. 


Table D: Case (4) Data. 


Value 

V 

of 

Node 

No. 

Max . 
Normal 
def . 

xl0~ 2 

(mm) 

Node 

No. 

Max . 
Tang . 
def . 

xlO -2 

(mm) 

Node Max. 

No . Compr . ’ 

stress 

(N/mm 2 ) 

Node 
No . 

Max . 

Shear 

stress 

(N/mm 2 ) 

0.02 


10 

0.4 1 

10 

0.386 

5 12.0 

5 

1 . 6 

0.2 


10 

0.440 

10 

0.385 

5 11.0 

5 

1.3 

0.4 


10 

0.440 

10 

0.384 

5 9.0 

5 

1.2 


The 

study 

of normal deflections 

(See Figure 

3.4) 

indicates 

that 

the 

normal 

deflections 

are high 

at the outer 

edge of the 

V- joint . 

On the side 

CD, their magnitude is higher 

and 

they rise 


more rapidly. Change in friction affects the normal deflections 
only over a small portion of the side CD. The deflections are 
found to reduce slightly with increase in . As the force is 
moved from location (1) to (4), the maximum value of normal 
deflection increases on both sides. 

The study of tangential deflection (See Figures 3.5 and 3.6) 
gives the information about the relative slipping between the two 
bodies. On side AB, it is the carriage which slips relative to 
the bed and that too in the downward direction, whereas on side 
CD, it is the bed which slips and it is in upward direction. The 
values of tangential deflection decrease marginally as p is 
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increased, almost uniformly over both the sides of the contact 
region . 

For case (1) (Figure 3.7(a)), it is seen that the maximum 
compressive stress is taken by node 6. Node 5 experiences less 
compressive stress. This difference gets reduced when the force is 
at location (2) (Figure 3.7(b)). At location (3) (Figure 3.7 
(c)), the peak compressive stresses at nodes 5 and 6 become almost 
equal whereas at location (4) (Figure 3.7 (d)), the maximum stress 
shifts to node 5. stresses at other nodes on side AB reduce 
progressively when the force moves from location (1) to (4). The 
shear stresses show a trend similar to that of the normal 
stresses. In each case it has been observed that with the 
increase in , the peak stress falls and the stress curves tend to 
flatten. 
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CHAPTER IV 


CONCLUSIONS AND SUGGESTIONS FOR FUTURE WORK 

4.1 CONCLUSIONS: 

The following conclusions can be drawn from the study of the 
planer bed-carriage contact problem. 

1. The method developed can be used to analyse any complex 
geometry, including unsymmetric geometry. As it uses 
substructuring technique, it can be used to solve very large 
problems. Very few iterations are required for convergence. The 
method can account for extremely low or high values of friction. 

2. Stresses in the carriage are very much less than the ultimate 
values but probably strength is not the criterion used for 
designing a carriage, it is the deflection, which affects the 
machining accuracy. 

3. Shear stresses which contribute to wear have been found to be 
larger at the bottom of the V joint. Differential use of costlier 
material or case carburising of such areas can extend the useful 
life of the machine. 

4. Stress concentration has been observed at the sharp corners 
of the joint. 

5. The maximum deflection observed seems to be within the 
acceptable machining accuracy level. However, stiffness members 
between the two V-slides of the bed may help in further improving 
the machining accuracy. 


42 



A.. 2 SUGGESTIONS FOR FUTURE WORK 


(i) The method can easily take care of variable coefficient of 
friction. However, experimental data is lacking for use in the 
study of contact problems. 

(ii) The method can be extended to tackle 3-D problems so that the 
problem of V-slide of a planer bed-carriage can be analysed more 
realistically. 

(ill) Multiple contact problems can be studied. 

( i v ) The m .thod presented can be modified for incremental loading 
by increasing the force from a low value to the required maximum 
in several steps. At each load step, it may be determined whether 
a node slips or adheres and appropriate contact conditions may 
thus be applied progressively. 

(v) Higher order elements can be used for better results. 
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